Two-dimensional electron system in high magnetic fields: Wigner crystal vs. 

composite-fermion liquid 
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The two dimensional system of electrons in a high magnetic field offers an opportunity to investigate 
a phase transition from a quantum liquid into a Wigner solid. Recent experiments have revealed an 
incipient composite fermion liquid in a parameter range where theory and many experiments had 
previously suggested the Wigner crystal phase, thus calling into question our current understanding. 
This Letter shows how very small quantitative corrections (< 1%) in the energy due to the weak 
interaction between composite fermions can cause a fundamental change in the nature of the ground 
state, thus providing insight into the puzzling experimental results. 
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When electrons confined in two dimensions are ex- 
posed to a strong magnetic field, they fall into the low- 
est quantum of their kinetic energy (Landau level), and 
the physics is governed by the interaction alone. With 
the kinetic energy thus quenched, one might have ex- 
pected electrons immediately to form a lattice, known as 
the Wigner crystal^, but experiments discovered the frac- 
tional quantum Hall effect^ instead, which is understood 
in terms of a correlated liquid of composite fermions'^"^ . 
The phase transition from this exotic quantum liquid 
into the electron solid anticipated at yet higher magnetic 
fields has been a subject of intense investigation. 

In a range of filling factor v, defined as the ratio of the 
number of particles to the Landau level (LL) degeneracy, 
interacting electrons in the lowest LL capture an even 
number of quantized vortices to transform into compos- 
ite fermions (CFs), because that is how they avoid one 
another most effectively. There are many types of com- 
posite fermions distinguished by their vorticity 2p, de- 
noted by 2pCFs. The effective filling factor of ^^CFs, v*, 
is given by the relation v = v* /{2pv* ± 1). (The CFs are 
often thought of as electrons carrying an even number of 
magnetic fiux quanta; see Fig. 1). It is often a good first 
approximation to treat the composite fermions as non- 
interacting. The experimental observation of fractional 
quantum Hall effect (FQHE) at v ^ n/{2pn ± 1) is well 
explained as the v* = n integral quantum Hall effect^ 
of ^PCFs, and extensive studies*"^ have shown that the 
model of noninteracting composite fermions is also quite 
accurate quantitatively. In particular, Laughlin's trial 
wave function® dX v — 1/m, where m = 2p -I- 1 is an 
odd integer, is interpreted as the ground state of non- 
interacting ^^'CFs ai V* = 1. Variational studies com- 
paring its energy with the energy of the Wigner crystal 
have found®"^'^ that the liquid wins for v = 1/Z and 1/5, 
but the WC has lower energy at v ~ 1/7 and smaller 
filling factors. Although the WC has so far eluded a 
direct observation, the insulating state observed in the 
region v < \/b has been interpreted as a Wigner crys- 
tal pinned by disorder^*"^^. However, transport^^ and 



photoluminescence^® experiments reported evidence for 
FQHE- like structure in the neighborhood of v = 1/7 
and V ~ 1/9 more than a decade ago. The recent 
observation^^, in exceptionally high quality samples, of 
a series of resistance minima at filling factors v — 1/7, 
2/11, 2/13, 3/17, 3/19, 1/9, 2/15 and 2/17 has given 
a particularly compelling indication for the existence 
of an incipient liquid of ^CFs and ®CFs in the range 
1/5 > 1/ > 1/9. 

While the approximation of neglecting the weak resid- 
ual interaction between the composite fermions is accu- 
rate, it is not exact. Because even small quantitative dis- 
crepancies can have substantial qualitative consequences 
near a phase boundary, we proceed below to calculate 
corrections due to the inter-CF interaction to determine 
how it modifies the previous results. 

The calculations will be performed in the spherical 
geometry, where N electrons are confined to the two- 
dimensional surface of a sphere, moving under the in- 
fiuence of a radial magnetic field produced by a mag- 
netic monopole of strength Q at the center. According 
to Dirac's quantization condition, Q can be either an in- 
teger or a half integer, and produces a total flux of 2Q(f>Q, 
where 0o = hc/e is called the flux quantum. Interact- 
ing electrons at Q map into weakly interacting ^^'CFs at 
Q* = Q — p{N — 1). We win investigate the n = 1 state 
of ^CFs, ®CFs, and ^°GFs, corresponding to FQHE at 
J/ = 1/7, 1/9, and 1/11, respectively. 

The system of noninteracting composite fermions is de- 
fined by analogy to noninteracting electrons at v* . In 
particular, the ground state of noninteracting composite 
fermions at i/* — n is obtained from the ground state 
of noninteracting electrons at v* = n, namely the state 
with n filled Landau levels, $„, as 



(1) 



where $i is the wave function of the lowest fully occupied 
Landau level and Plll is the lowest LL projection oper- 
ator. The factor <i>J^ attaches 2p vortices to each electron 
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to convert it into a composite fermion, so the wave func- 
tion on the right hand side may be viewed as n filled levels 
of composite fermions. For n = 1, the above wave func- 
tion, shown pictorially in Fig. la, reduces to Laughlin's 
wave function. The residual interaction between com- 
posite fermions modifies 5'^°-' through hybridization with 
higher CF levels. We consider 

^-0) = P^^^$2p^o) (2) 

where denote states obtained from by exciting 
j electrons from the topmost occupied LL to the low- 
est unoccupied LL, with a labeling distinct particle-hole 
configurations, ^'a ' thus contain j particle-hole pairs of 
composite fermions (Fig. 1). An improved representa- 
tion for the ground state, x^''^\ which involves the effect 
of CF level mixing caused by the inter-CF interaction, is 
obtained by diagonalizing the Coulomb Hamiltonian in 
the basis given by [^(o), {*^^'}, ...{4'^"'^}]. The 

corresponding energy per particle, E^'^\ gives a strict 
variational upper bound for the true ground state en- 
ergy. Successively higher values of J produce lower en- 
ergies, and convergence is achieved when E^'^^ becomes 
more or less independent of J. 

With explicit wave functions in hand, the computation 
is conceptually straightforward, but technically challeng- 
ing due to the complex form of wave functions. Several 
complications are handled as follows. The lowest LL pro- 
jection is evaluated as discussed in the literature^*^, by re- 
placing the single particle eigenstates in by appropri- 
ate wave functions for composite fermions. In the spheri- 
cal geometry, the liquid state, which possesses rotational 
symmetry, has total orbital angular momentum L = 0, 
as appropriate for a filled angular momentum shell. The 
Coulomb interaction does not mix states with different 
L, so it is convenient to work within the L = sector. 
We generate the many body states ' in this sector by 
diagonalizing the Coulomb interaction at v* = n in the 
subspace containing j particle-hole pairs. The resulting 
eigenstates are a linear superposition of Slater determi- 
nant basis states, from which basis functions at v are 
obtained according to Eq. (2). 

The angular momentum of each successive LL increases 
by one unit in the spherical geometry. For the states with 
a single particle hole pair, the angular momenta of the 
particle and the hole differ by one unit, and as a result fail 
to produce a state with total angular momentum L = 0. 
It is therefore necessary to consider J > 2. We consider 
in this work up to three particle hole pairs (J = 3); for 
still higher values of J, the basis becomes too large to 
produce numerically stable energies. Our inability to go 
beyond J = 3 imposes a limitation on the system size for 
which we can obtain reliable results. 

Another major obstacle is that the basis functions at 
Q generated from the wave functions at Q* are not or- 
thogonal. The standard Gram-Schmidt procedure will 
be employed to construct an orthogonal basis. The inner 



products needed for this purpose, as well as the matrix 
elements of the Coulomb Hamiltonian needed for diago- 
nalization are computed by Monte Carlo. Approximately 
10 Monte Carlo runs were performed for each energy, 
with 0.5-1.0 xlO^ iterations in each run. The calculation 
was performed on a PC cluster, and took the equivalent 
of > 100 days of CPU time on a single node (dual 1 
GHz Intel PHI). The statistical uncertainty in each ma- 
trix element translates into an uncertainty in the energy, 
which we estimate by calculating it in ~ 10 different runs. 
Satisfactory accuracy is obtained with the number of it- 
erations performed. More details on the method can be 
found elsewhere^ ^. All energies quoted below include in- 
teraction with a uniform, positively charged background, 
and have been multiplied by a factor ^2Qu/N to correct 
for the finite size deviation in the density. A strictly two 
dimensional system is assumed and a mixing between the 
LLs of electrons is neglected. 

There are substantial fluctuations in energy for N < 
6, which we believe arise from commcnsurability effects. 
For = 4 and 6, polyhedral arrangements are possible 
on the sphere (tetrahedran and octahedran) producing 
anomalously low energies. For A'^ = 3 there is only a 
single L = state at = 1/m, indicating that the system 
is too small to be meaningful. In the following, we only 
consider N > 7, where we have found that the ground 
state energy behaves more or less smoothly with 1/N. 
(Polyhedra are possible also for certain larger values of 
A'', but there the quantitative effect is small.) 

We begin by testing the accuracy of the method. Table 
1 shows and E'-^^ for several A^ at ly = 1/5, along 
with exact energies obtained by numerical diagonaliza- 
tion. E^"^^ has a remarkable accuracy of 0.01% or better, 
which is a significant improvement over which is off 
by up to 0.2%. Fig. (2) shows the energies E("\ E^'^\ 
and as a function of N-^ for V = l/l, 1/9, and 
1/11. For N = 7, 8, and 9, convergence is manifest. For 
N = 10, E^'^^ is not yet J independent, but we expect 
to be reasonably accurate. Larger values of J must 
be considered with increasing N in order to obtain mean- 
ingful estimates; for a fixed J, the difference between the 
energies (per particle) E'-^^ and E'-'''^ is 0{N-^), which 
vanishes in the thermodynamic limit. With J < 3, we are 
able to obtain reliable estimates for the energy for only 
up to A" = 10. (It is noted that the dimension of the 
full Fock space for 10 particles at = 1/11 is 1.6 x 10^^, 
rendering exact diagonalization impracticable.) 

We use a linear extrapolation in A^~^ to obtain the 
thermodynamic limit (solid line), as appropriate for 
a compact geometry without boundaries (boundaries 
would contribute a term ^ N^^^^ in the finite size cor- 
rection for the energy per particle). To get a feel for 
the accuracy of such an extrapolation with four points, 
we note that for E^^\ which serves as a useful bench- 
mark, the linear extrapolation using the first four points 
(dotted line) differs from the more accurate extrapolation 
(dashed line) by approximately 0.04%, 0.1%, and 0.15% 
for v = 1/7, 1/9, and 1/11, respectively. 
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E<-"\ the energy of Laughlin's wave function, exceeds 
the WC energy for ly = 1/7, 1/9, and 1/11, as noted 
in earher studies^' However, the best variational es- 
timate, is ~ 0.8% lower than for the filling 
factors considered, which is not a large reduction in ab- 
solute terms, but enough to make a qualitative change 
in the outcome. The principal conclusion of this work 
is that, for the model considered, the ground state is a 
liquid at = 1/7 and very likely also at ly = 1/9. At 
u = 1/11 the relatively large uncertainties in various en- 
ergies prevent us from making a definite statement, but 
a liquid ground state is a distinct possibility even here. 
It is striking that composite fermions with vorticity as 
high as eight or perhaps even ten can occur in nature. 

The relatively small overlaps between \E'(°) and x^^^ 
(Fig. 2) are indicative of a substantial amount of mixing 
between CF levels. Nevertheless, there is little doubt that 
the true ground state can be obtained from the noninter- 
acting state ^^^^ without crossing a phase boundary (or 
without gap closure), and hence the physics of the true 
liquid state can qualitatively be described in terms of non- 
interacting composite fermions. Only detailed and accu- 
rate quantitative calculations can tell, however, whether 
the liquid or the solid is the ground state. 

The local resistance minima observed at fractions of 
the form i/ = n/ (6n ± 1) and v = nl (8n ± 1) ride a back- 
ground that rises rapidly with decreasing temperature^^. 
Pan et al}^ suggested that the CF liquid is obtained 
at temperatures exceeding the WC melting temperature. 
Our calculations put forth the alternative possibility that 
the true gromid state at these fractions, in the absence 
of disorder, is the CF liquid. The insulating behavior 
would then be the result of an absence of percolation 
of the liquid, due to domains of Wigner crystal induced 
by disorder^^. A future observation of a well developed, 
dissipationless FQHE state at these fractions would con- 
firm this scenario. (It may be recalled that the FQHE 
at = 1/5 was also initially observed on a rising back- 
ground.) Interesting questions also arise regarding the 
nature of the state away from these special filling factors. 
It is possible that there is a series of zero-temperature 
quantum phase transitions between the CF liquid and 
the WC as a function of i^, as is believed to be the case^** 
in the vicinity oi v = \ jh. 

To summarize, the interaction between composite 
fermions is weak, causing only a small change in the en- 
ergy (< 1%), but it can be of crucial importance in sta- 
bilizing the CF liquid in parts of the phase diagram ear- 
lier thought to belong to the Wigner crystal. We thank 
Vito Scarola for help and useful discussions. This work 
was supported in part by the National Science Founda- 
tion under grants no. DGE-9987589 (IGERT) and DMR- 
0240458. We are grateful to the High Performance Com- 
puting (HPC) Group led by V. Agarwala, J. Holmes, and 
J. Nucciarone, at the Penn State University ASET (Aca- 
demic Services and Emerging Technologies) for assistance 
and computing time with the LION-XE cluster. 
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TABLE I. The ground state energy per particle for N elec- 
trons at V — 1/5. The Laughlin energy E^°^ (also the energy 
of the noninteracting CF ground state), and E^'^\ which in- 
cludes the effect of inter-CF interaction at the lowest level of 
approximation, are compared with the exact energy obtained 
from numerical diag onaliza tion. The energies are quoted in 
units of e'^ /el {I = ^ Tic/eB is the magnetic length and e is the 
dielectric constant of the host semiconductor) , and the statis- 
tical uncertainty from the Monte Carlo simulation is shown 
in parentheses. For three particles, E^'^^ and E^'^^ are exact 
for the trivial reason that there is only one uniform state. 
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FIG. 1. Examples of configurations of ^CFs at v* — 1 
[v — 1/3). Panel (a) shows pictorially the ground state 
of noninteracting composite fermions, while panel (b) shows 
two ^CF particle-hole pairs. The analogous configurations at 
v = 1/7, 1/9, and 1/11 will contain ^CFs, ®CFs, and ^°CFs, 
respectively. 
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FIG. 2. The ground state energy per particle as a function 
of the number of composite fermions, A'^, for v — 1/7,1/9, 
and 1/11. E^'^^ is the energy of the noninteracting CF sys- 
tem, while E^^^ and E^^^ are successively better approxima- 
tions including interaction effects. The cross on the y-axis 
shows the energy of the Wigner crystal, taken from Ref. 9, 
which has an uncertainty of ~ 0.0005 in the chosen units. 
The error bars denote the statistical uncertainty in the Monte 
Carlo and the solid straight lines are the best chi-square 
fits for E^^\ The number near each solid diamond equals 

I < x'^'l*'") > |V(< x<'^!x<'' X ^'^'l*'"' >). 
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